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In this paper, the holographic dark energy in Brans-Dicke theory is confronted by cosmic observa- 
tions from SN la, BAO and CMB shift parameter. The best fit parameters are found in la region: 
n h0 = 0.683lg;gi, c = 0.605±g;ig| and a = 0.006621^67 (equivalently u = 905.690l^;t?? which 
is less the solar system bound and consistent with other constraint results) . With these best fit val- 
ues of the parameters, it is found the universe is undergoing accelerated expansion, and the current 
value of equation of state of holographic dark energy Who = — 1.246 + g '{44 which is phantom like in 
Brans-Dicke theory. The evolution of effective Newton's constant is also explored. 
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1 1' I. INTRODUCTION 

o\ 

The observation of the Supernovae of type la [J, 0] provides the evidence that the universe is undergoing accelerated 
expansion. Jointing the observations from Cosmic Background Radiation HH and SDSS one concludes that 

the universe at present is dominated by 70% exotic component, dubbed dark energy, which has negative pressure and 
1 ■ push the universe to accelerated expansion. Of course, the accelerated expansion can attribute to the cosmological 
constant naturally. However, it suffers the so-called fine tuning and cosmic coincidence problem. To avoid these 
problem, dynamic dark energy models are considered, such as quintessence B , H, [ 9lll0ll ■ p htantom pd| . quintom [ijj 
and holographic dark energy [l3|, [3] etc. For recent reviews, please see [Til [la. 17, 18[ [ljl [2(j ■ In particular, the 
—t, model, named holographic dark energy, is constructed by considering the holographic principle and some features 
of quantum gravity theory. According to the holographic principle, the number of degrees of freedom in a bounded 
system should be finite and has relations with the area of its boundary. By applying the principle to cosmology, one 
can obtain the upper bound of the entropy contained in the universe. For a system with size L and UV cut-off A 
without decaying into a black hole, it is required that the total energy in a region of size L should not exceed the 
mass of a black hole of the same size, thus L 3 p^ < LM^. The largest L allowed is the one saturating this inequality, 
thus pa = 3c 2 A/p Z L -2 , where c is a numerical constant and M p i is the reduced Planck Mass M~f = 8nG. It just 

■ means a duality between UV cut-off and IR cut-off. The UV cut-off is related to the vacuum energy, and IR cut-off is 
related to the large scale of the universe, for example Hubble horizon, event horizon or particle horizon as discussed 

■ by [H, [l4| . In the paper 14 1 , the author takes the future event horizon 
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as the IR cut-off L. This horizon is the boundary of the volume a fixed observer may eventually observe. One is to 
formulate a theory regarding a fixed observer within this horizon. As pointed out in (l4| . it can reveal the dynamic 
nature of the vacuum energy and provide a solution to the fine tuning and cosmic coincidence problem. In this model, 
the value of parameter c determines the property of holographic dark energy. When c > 1, c = 1 and c < 1, the 
holographic dark energy behaviors like quintessence, cosmological constant and phantom respectively. 

On the other hand, Brans-Dicke theory [2ll as a natural extension of Einstein's general theory of relativity can pass 
the exp erimental tests from the solar system [22( and provide explanation to the accelerated expansion of the universe 
[25L 12a , I27I ] . In Brans-Dicke theory, the gravitational constant is replaced with a inverse of time dependent scalar field, 
i.e. 8nG = tp^y, which couples to gravity with a coupling parameter u>. The holographic dark energy model in the 
framework of Brans-Dicke theory which has already been considered by many authors [2l|, [H, [H, H3, HH, HH . In our 
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previous paper [2l| . the properties of holographic dark energy in Brans-Dicke theory was discussed by giving some 
characteristic values of the parameters, where the values of the parameters were given by taking the corresponding 
values obtained from constraint by cosmic observations in Einstein theory. However, in Brans-Dicke theory, the values 
of parameters would be different from that in Einstein theory. After all, the Newton's constant G evolves with time in 
Brans-Dicke theory. So, there would be some differences. In fact, in our previous paper [2l[, the value of \u>\ > 40000 
is taken for granted. It would be dangerous because of the possibility of small value in large scale, say in cosmological 
scale, reported by the authors [34| . So, the holographic dark energy in Brans-Dicke theory must be tested by cosmic 
observations. This is the main task of this work. In this paper, the cosmic observations from SN la, BAO and CMB 
shift parameter will be used as cosmic constraints, for details please see the following sections. When using these 
observational data set, one has to notice the evolution of Newton's constant G. In fact, the SN la as a useful cosmic 
constraint has been considered in [34|, [H, [H, H3] ■ 

This paper is structured as follows. In Section [III we give a brief review of holographic dark energy in Brans-Dicke 
theory. A brief description of comic observations and methods used in our paper is listed in Section IIIIl Section IIVI 
is the results and discussion. 



II. HOLOGRAPHIC DARK ENERGY IN BRANS-DICKE THEORY 



Here, we just give a brief review of holographic dark energy in Brans-Dicke theory, for the details please see [3c 
The holographic dark energy in Brans-Dicke theory takes the form 



Ph = 3c 2 $(i)L- 2 , 



(2) 



where 



is a reverse of time variable Newton's constant. In a spatially flat FRW cosmology filled with 
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dark matter and holographic dark energy, the gravitational equations can be written as 



3$ 



2 $ LU& 



Ph, 



2- + H 2 + -— + 2H- + - = 



(3) 
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where H = ^ is the Hubble parameter, p m is dark matter energy density, ph is the holographic dark energy density 
and pu is the pressure of holographic dark energy. The scalar field evolution equation is 



$ + 3H<I> = 



Pm + Ph~ 3Ph 



(5) 
(6) 



2uj + 3 

Considering the dark matter energy conservation equation 

Pm + 3Hp m = 0, 

and jointing it with Eq. Q, Eq. ([4]) and Eq. (O, one obtains the holographic dark energy conservation equation 

Ph + 3H(p h +Ph) = 0. (7) 
Here, we have considered non-interacting cases. The Fricdmann equation (J3J) is 



H 2 = Pn 



£lL- H- + - — 



With the assumption $/<E>o = ( a /oo) Q ! the Eq. ([8]) is rewritten as 

2 



H 1 



(6 + 6a - wa 2 )$ 



{pr, 



Ph 



(8) 



(9) 



It is easy to find out that, in the limit case a — * 0, the standard cosmology is recovered. To make the Friedmann 
equation ((9]) to have physical meanings, i.e. to make (6 -I- 6a — wa 2 ) > 0, one has the following constraints on the 
values of a 



<a< for w > 



a < 



or a > 
ft, 



3 +v /9+6^ for _3/ 2 < u < 0, 

for u < -3/2. 



(10) 
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However, the solar system experiments predict the value of a; is |u>| > 40000 [22j. However, the value of parameter 
u) = —3/2 is a boundary of ghost [23[. So, in this paper, when considering these constraints, the second line of Eq. 
(fTOf will be omitted and u> > will be consider in this paper. In fact, authors [24[ have used the cosmic observations 
to constrain the parameter ui. In [24| . the authors found that to can be smaller than 40000 in cosmological scale, say 
w ~ 1000. 

When the event horizon R e h 

f°° dt f°° da 

R ^- a J t W) =a Ja 

is taken as the IR cut-off. The holographic dark energy is 



3c 2 $ 
R lh 

And, the Friedmann Eq. ^ is rewritten as 



Ph = -s-. (12) 



h 2 = Hfo m0 (^) (3+a) + il h H 2 



= H 2 n mQ a-^ +fl h H 2 , (13) 

of dark matter is fl m 
holographic dark energy Qh is the solution of differential equation 



where the dimensionless energy density of of dark matter is fi m o = (s+sa-ua 1 ) <j^jfs = 87r 3 ^r Q an d the one of 



where ' denotes the derivative with respect to x = In a. This equation describes the evolution of dimensionless energy 
density of dark energy. Comparing the definition of Q m0 with that of the standard cosmological model, one easily has 
the relation toa = 6. With the relation ao/a = 1 + z, the Friedmann equation ()13|) can be rewritten as 

H =H ° i-n h =H ° rr^ • ( 15 ) 

From the conservation equation of the holographic dark energy ([7]), on has the equation of state (EoS) of holographic 
dark energy 

where Wh = Ph/ Ph- From the above equation, one finds the EoS of holographic dark energy is in the range of 

-^(l + a+^j <w h <-^(l + a), (17) 

when one considers the holographic dark energy density ratio < flh < 1. Also, by using the Eq. Q and the 
assumption §/&o — (a/a>o) a , one obtains the deceleration parameter as follows 

a 1 a 6w h n h rie , 

q = -^p = 2 +a+ 8T^ + ^T^- (18) 

It is clear that the 'Standard' holographic dark energy will be recovered in the limit a — * 0. In Brans-Dicke theory 
case of holographic dark energy, the properties of the holographic dark energy are determined by the best fit values 
of parameters c and a which would be obtained by confronting with cosmic observational data. It can be easily seen 
that the holographic dark energy can be quintessence, phantom and quitom as that in the Standard case. But, all 
these properties must be determined by cosmic observations. 



III. COSMIC OBSERVATIONAL CONSTRAINTS 



In this section, cosmic observations and methods used in this paper are described. 
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A. SN la 



We constrain the parameters with the Supernovae Cosmology Project (SCP) Union sample including 307 SN la 
[HI, which distributed over the redshift interval 0.015 < z < 1.551. Constraints from SN la can be obtained by fitting 
the distance modulus fj,(z) [U, l35l. 1361 |37| 



15 G 

H th (z) = 5 log 10 (D L (z)) + — log 10 + no, 



(19) 



where, G is the current value of effective Newton's constant G e ff, Dl(z) is the Hubble free luminosity distance 
Hq<1l{z) / c and 



d L (z) = c(l + z) 



G dz' 



(20) 



/„ \G{z')H(z') 

Ho = 42.38- 51og 10 ft, (21) 

where H is the Hubble constant which is denoted in a re-normalized quantity h defined as H — 100ft, km s _1 Mpc _1 . 
The observed distance moduli fJ-obs(zi) of SN la at Zi is 

Hobs(zi) = m obs {zi) - M, (22) 

where M is their absolute magnitudes. 

For SN la dataset, the best fit values of parameters in a model can be determined by the likelihood analysis is based 
on the calculation of 

X {Ps,m a ) = 



SNIa 



E 

SNIa 



[5 log 10 (D L (p s ,Zi)) - m obs (zi)+m y 



(23) 



where mo = (J-o + M is a nuisance parameter (containing the absolute magnitude and Ho) that we analytically 
marginalize over [39l ]. 



X 2 (p s ) - -2 In 



to obtain 



where 



exp 



B 



1 



X = A - — + In — 



X 2 (p s ,rn ) 



G 



drrio 



C 



2tt 



SNIa 

B = E 

SNIa 



[51og 10 (Di(p a ,^)) - m obs(zi)Y 



5 log 10 (D L (j>s ,Zi) - m obs (zi 



SNIa 1 



(24) 

(25) 

(26) 
(27) 
(28) 



The Eq. 



(29) 



has a minimum at the nuisance parameter value mo = B/C. Sometimes, the expression 

X 2 S Nia(Ps,B/C)=A-{B 2 /C) 

is used instead of Eq. (|25[) to perform the likelihood analysis. They are equivalent, when the prior for mo is flat, as 
is implied in (|24|) , and the errors Oi are model independent, what also is the case here. Obviously, from the value 
mo = B/C, one can obtain the best-fit value of h when M is known. 

To determine the best fit values of parameters for each model, we minimize \ 2 (j> s ,B jC) which is equivalent to 
maximizing the likelihood 

C( Ps ) oc e -x-( Ps ,B/c)/2^ (30) 
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BAO 



The BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure 
the distance-redshift relation at z = 0.2. BAO in the clustering of the SDSS luminous red galaxies measure the 
distance-redshift relation at z = 0.35. The observed scale of the BAO calculated from these samples and from the 
combined sample are jointly analyzed using estimates of the correlated errors, to constrain the form of the distance 
measure D v (z) EJ ED, El 



D v (z) 



(l + z) 2 D 2 A (z)- 



1/3 



(31) 



H(z) 

where Da(z) is the proper (not comoving) angular diameter distance which has the following relation with d^iz) 

dm - W) 

Matching the BAO to have the same measured scale at all redshifts then gives El 

Dy(0.35)/£V(0.2) = 1.812 ±0.060. (33) 

Then, the X%Ao(Ps) is given as 

2 , x [£>v(0.35)/Dy(0.2) - 1.812] 2 
Xbao(Ps) = ■ (34) 

C. CMB shift Parameter R 

The CMB shift parameter R is given by El 



R{z*) = ^n m H$(l + z*)D A (z*)/c (35) 



which is related to the second distance ratio D A (z*)H(z*)/c by a factor y/1 + z*. Here the redshift z* (the decoupling 
epoch of photons) is obtained by using the fitting function [44j | 

z, = 1048 [l + 0.00124(fVi 2 )- - 738 ] [1 + gi(n m h 2 ) 92 ] , (36) 

where the functions g\ and gi are given as 

9l = 0.0783(n b h 2 )-°- 23S (l + 39.5(n b h 2 )°- 763 y\ (37) 
.92 = 0.560(1 + 21. l(Q b h 2 ) 1 - 81 y 1 . (38) 



The 5-year WMAP data of R(z*) = 1.710 ± 0.019 [45J will be used as constraint from CMB, then the Xcmb(Ps) is 
given as 

2 , (R(zJ - 1.710) 2 
X 2 cmb(Ps) = Q0192 (39) 

IV. RESULTS AND DISCUSSION 

For Gaussian distributed measurements, the likelihood function L oc e~ x / 2 , where x 2 is 

X 2 = XsNIa + Xb.AO + XCMB, (40) 

where XsNIa ^ s gi yen m Eq. (|29]) . Xbao * s gi yen m Eq. p4|) . Xcmb ^ s given in Eq. ([39]) . In this paper, the 
central values of Q b h 2 = 0.02265 ± 0.00059, Q m h 2 = 0.1369 ± 0.0037 from 5-year WMAP results [H] and H = 
72 ± 8kms _1 Mpc _1 are adopted. After calculation, the results arc listed in Tab. [B The contours of flho — c and 
&ho ~ c with la, 2a confidence levels are plotted in Fig. [1] 
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Datasets 




^feo(lc) 


C(lff) 


a(lff) 


uj(la) 


SN+BAO+CMB 


317.156 


0.6831°;™ 


0.6051°,;^ 


00 6 62+ 00477 

u - UUDUZ -0. 00467 


905.6901SI3?; 



TABLE I: The minimum values of x 2 and best fit values of the parameters. 
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FIG. 1: The contour plot of Qh.o — c and flho — ot with la, 2a confidence levels. 



Current constraints 46] on the variation of Newton's constant imply 

G eff 



G 



eff 



< 10- n yr-\ 



in our case, which corresponds to 



It implies 



aH < 10~ n yr~ l . 



a < — x 10 _11 w _1 - 
~ H y 



(41) 



(42) 



(43) 



Considering the current value of Hubble constant h = 0.721°;^, one obtains the bounds on a, when the central value 
is taken 



a < 0.135807. 

It is clear that the best fit value of parameter a is under the bound and consistent. 
The evolution of the effective Newton's constant is written as follows 



G 



G 



= (1 + *Y 



(44) 



(45) 



under the assumption $/$o = ( a / a o) a - With the best fit value of the parameter a — 0.00662, the evolution of 
effective Newton's constant with redshift z is plotted in Fig. [2] with the best fit parameter. 



V. CONCLUSIONS 



In this paper, the holographic dark energy in Brans-Dicke theory is constrained by cosmic observations which 
include the data sets from SN la, BAO and CMB shift parameter. In ler region, the best values of the parameters 
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,3 O 1.Q4 - 




FIG. 2: The evolution of effective Newton's constant with redshift z, with la regions. 



are: Q h0 = 0.683±g:g||, c = 0.605t£^ and a = 0.00662±g;[$^. Equivalently, with the relation cua = 6 holds, the 
best fit value of u> in la region is to = which is smaller than the value from the solar system bound, 

but consistent with other reports in cosmological scale 34J. With these best fit values of parameters, it is found the 
universe is undergoing accelerated expansion currently, and the current value of equation of state of holographic dark 
energy Who = — 1. 246^' 144 which is phantom like in Brans-Dicke theory. The evolution of effective Newton's constant 
is explored. Its best fit value is consistent with the bound [461 ] . 
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